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Abstract 

An asymmetric covering T>(n, R) is a collection of special subsets S of an n-set such that every 
subset T of the re-set is contained in at least one special S with \S\ — \T\ < R. In this paper 
we compute the smallest size of any T>(n, 1) for re < 8. We also investigate "continuous" and 
"banded" versions of the problem. The latter involves the classical covering numbers C(n, k,k — l), 
U ; and we determine the following new values: C(10,5,4) = 51, C(ll,7,6,) = 84, C(12,8,7) = 126, 

C(13, 9,8) = 185 and C(14, 10, 9) = 259. We also find the number of nonisomorphic minimal 
covering designs in several cases. 

1. Introduction 

£> \ Let D(n,R) denote the smallest size of any asymmetric covering^ T>(n, R). Prompted by applica- 

tions to the manufacture of semiconductor wafers, Cooper, Ellis and Kahng |7|] have investigated 
£f~^ , the asymptotic behavior of D(n,R) for fixed R as re — > oo. In Section 2 of the present paper we 

IQ | show that the values of D(n, 1) for n < 8 are as shown in Table [|{] 

Our method of attack is to formulate D(n, 1) as the solution to a {0, l}-integer programming 



o 

■ problem. If instead we allow the variables to take any real values in the range [0, 1] then the linear 

program can be solved exactly (Section 3). Of course this provides a lower bound E(n) to D(n, 1). 

An upper bound can be obtained by restricting to asymmetric coverings with a certain banded 
structure defined in Section 4. Corollary [2] shows that the solution C(n) to the banded version of 
the problem is given by 

> ■ [n/2] 

j><J : C{n) = C(n + l,n + l-2i,n-2i) , (1) 

H ' i=0 

where as usual C(v,k,t) denotes the smallest size of any covering design C(v,k,t); that is, any 
collection of special fe-subsets S of a v-set such that any t-subset T is contained in at least one S. 

Although there have been a large number of papers written about covering designs (for recent 
work see f|, 1 1C ] , [11], |^5| , fill , [17], [19]), not many exact values are known. In Section 5 we 



determine several new values of C(n, k,k — l). Call a covering design C(n, k, t) optimal if it contains 
the smallest number C(n,k,t) of subsets, and minimal if it is no longer a covering if any subset 
is omitted. In Section 5 we also determine the number of nonisomorphic optimal covering designs 
C(n,k,k — 1) for n up through 10, as well as the number of minimal covering designs for n up 
through 7. See Tables [|, ||| in Section 5. 

1 Present address: Center for Communications Research, Princeton, NJ 08540 

2 These are called directed coverings in However, that term has already been used in the literature with a 
different meaning (cf. pB|). 

3 This is sequence A66000 in Jl8[. If any further values are computed they will be recorded there. 
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Using these results we obtain the values of C(n) for n < 11 shown in Table [T]. 

Some of the literature on covering designs works with the complements of the special sets, in 
which case this is called the Turdn design problem (cf. Q), and the Turdn number T(v,k,t) is 
equal to C(v , v — t, v — k). Of course our results also provide new values for certain Turan numbers. 

Table 1: Size D(n, 1) of smallest asymmetric 1-covering of an n-set, together with values of the 
continuous and banded solutions E(n) and C(n). 
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Notation 

Let denote the set of binary vectors of length n. We represent subsets of an n-set by their 
indicator vectors in F^, and then an asymmetric covering T>(n,R) can be thought of as a binary 
code called an "asymmetric covering code" . As usual weight (denoted wt) and distance (dist) refer 
to Hamming weight and Hamming distance. The co-weight of u & Fj is n — wt(u). Two codes, 
coverings or designs are isomorphic if they differ just by a permutation of the coordinates. 

2. Values of D(n, 1) for n < 8. 

Let x u , u 6 Fg , be real {0, l}-valued variables. Then D(n, 1) is equal to the minimal value of 
subject to the constraints 

^2 x u + x v > 1, for all v G F£ , (3) 

where u D v indicates that u covers v and dist(u,v) = 1. If w = 11 ... 1 then necessarily x u = 1. 
The corresponding asymmetric covering code consists of the vectors u 6 FJ, for which x u = 1. For 
example D(3, 1) = 3, and the code (which is unique up to permutation of the coordinates) is 

{111,110,001} . (4) 
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Every binary vector of length 3 is either in this code or is contained in a codeword at distance 1 
below it. 

We call the above minimization problem the exact integer programming (or IP) problem. If we 
relax the constraints and allow the x u to take any real values in the range [0, 1] we get a continuous 
linear programming (or LP) problem, whose solution we denote by E{n). 

Theorem 1. The values of D(n, 1) for n < 8 are as shown in Table 0. 

Proof. We attacked the IP problem using CPLEX |8| with AMPL H as a convenient interface. 
CPLEX uses a branch and bound strategy for such problems. We regard solutions obtained in 
this way as perfectly rigorous, since the computations could in principle be replaced by extremely 
tedious hand calculations. 

For n < 7 CPLEX was able to find solutions directly, without any additional assumptions being 
added. Explicit solutions are described in Section 4. 

For n = 8 we must show that D(8, 1) = 58. A solution of size 58 found by CPLEX is given in 
Table |. (Each vector is represented by two hexadecimal characters. This covering has no apparent 
structure-in particular it has trivial automorphism group) To show that 57 is impossible we argue 
as follows. If we add the extra assumption that there are at most eight codewords of weight 5 
to the continuous LP problem, the solution is at least 63. Therefore there must be at least nine 
codewords of weight 5. From the tables of constant weight codes g it follows that there must 
be two codewords of weight 5 and distance exactly two apart. Without loss of generality we can 
assume that u\ = 11111000 and U2 = 11110100 are in the code. 

Suppose the code contains a vector U3 of weight 3 with dist{u\,u^) = 6 and dist{u2,u^) = 4, 
say U3 = 00010101. Then CPLEX finds that the minimal solution to the IP problem is 58. On the 
other hand if no such vector 113 is present (this rules out 16 vectors of weight 3) no feasible solution 
to the IP problem of size < 57 exists. Hence D{8, 1) = 58. ■ 

The total computing time for these calculations was less than 48 hours. 

Table 2: A minimal asymmetric covering T>(8, 1) of length 8 containing 58 sets (represented in 
hexadecimal) . 

01 07 OA 11 IE 28 2D 33 34 37 3B 4B 4C 52 55 57 5D 61 66 6E 
6F 73 75 78 7E 7F 84 89 8F 96 98 99 9F A2 A5 AA B3 BB BC BD 
CO C3 CC D5 DA DB DD E6 E7 E9 EE EF F0 F6 F7 F9 FE FF 



3. Solution to the continuous linear programming problem 

Theorem 2. The optimal solution to the continuous LP problem of choosing < x u < 1 foru € 
so as to minimize $y subject to ^) is given by 

E(n) = {-l) n n\{R n (2) - 1^(1)^(1)} (5) 

where 

" (_ T \k 



k\ 

k=0 



is the degree n partial sum of e x . 
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Proof. We may assume that x u depends only on the weight of u. (For let y u denote the average 
value of x v over all v with wt(v) = wt(u). Then by averaging ||) we see that the y u satisfy the 
same constraints as the x u , and 

^ ] Uu = ^ ^ %u ■ 

So a symmetrized solution is just as good as a general solution.) 

The "weight enumerator" of a symmetrized solution x u is defined by A w = Yl w t{u)=w x « f° r 
w = 0, . . . , n. The quotes are needed because the x u are in general not integers. Let A(z) = 

The covering condition (|3[) reads 

wA w + A w -i>[ U ), w = l,...,n, 
\w-l J 

or in other words 

A(z)+A'{z) > (z + l) n . (7) 

We wish to choose Ao,...,A n > so as to minimize A(l) subject to (0). The dual problem 
(compare [[U], Chapt. 17]) is to choose Bq, . . . , B n > so as to maximize B(l) subject to 

B(z)+B'(z)<(z + l) n (8) 

where B{z) = Y^=q B w z w . We claim that 

B(z) = {-l) n n\{R n (z + 1) - R n ^(l)R n (z)} (9) 

is a feasible solution to the dual problem. In fact it is straightforward to verify that B w > for all 
w and 

B(z) + B'(z) = (z + l) n - R n ^(l)z n . 

Since i?„_i(l) > 0, (|) holds. 
Therefore 

B(l) = (-l) n n\{R n (2) - Rn^Rn^l)} 

is an upper bound to the optimal solution to the primal problem. 
On the other hand 

A(z) = {-l) n n\{R n {z + 1) - Rn(l)Rn-l (*)} 

satisfies 

A(z) + A\z) = (z + l) n + R n {l)nz n - 1 

and is easily checked to be a feasible solution to the primal problem. Since A(l) = B(l), this must 
be the optimal solution to both problems. ■ 



Corollary 1. As n — ► oo, 



S(n) ~ 2"+! ( i - A + o f J_ ) ) . 



n v? \ n 3 



We omit the routine derivation of this from (||). 

The first few values of E(n) are shown in Table [l]. 
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4. Banded solutions 



Let T> be an asymmetric covering T)(n, 1). We call V banded if every vector v € F?> with odd 
co- weight is covered by a vector u £ D of weight one higher. 

For example (f|) is banded, since the vector 000 is covered by 001 and the vectors 011, 101, 110 
are all covered by 111. 



Theorem 3. // a code C C ¥^ +1 is a union of covering designs, 

[n/2] 

C = (J C(n + 1, n + 1 - 2i, n - 2t) , (11) 

i=0 

i/ien deleting^ any one coordinate from all the vectors of C yields a banded asymmetric covering 
T>(n, 1). Conversely, let T> be a banded asymmetric covering T>(n, 1). // we append a or 1 to every 
vector of T> in such a way that all co-weights become even, the result is a union of covering designs 



of the form ( 11 ). 



Proof. Suppose C has the structure shown in (11) and let T> be obtained by deleting one coordi 



nate, which for concreteness we suppose is the last coordinate. We must show that T> is a banded 
asymmetric covering. Let v € F?J have weight w. If the co- weight n — w is even, say 2i, then v* = vO 
must be covered by some vector u* = u5, 5 = or 1, in the covering design C(n + 1, n+1 — 2i, n — 2i), 
and then u € D covers v. On the other hand if n — w is odd, say 2i — 1, then v* = vl must be 
covered by some u* = ul € C(n + 1, n + 1 — 2i,n — 2i), and again u covers v. The converse is 
established by similar arguments. 

Since the covering number C(n, k, v) is by definition the size of the smallest C(n, k, v), we have: 

Corollary 2. The size of the smallest banded asymmetric covering T>(n, 1) is given by 

[n/2] 

C(n) = G{n + 1, n + 1 - 2i, n - 2%) . 

i=0 

Using the known values of C(n,k,t) and the new values to be established in the next section 
(see Table [|) we can determine C(n) exactly for n < 11. These values are given in Table and 
show that for n < 7, banded asymmetric coverings are as good as any asymmetric coverings. 

A more detailed investigation provides further information: 

Theorem 4. For lengths n = 1,2,3,5,6 and 7 an optimal asymmetric covering is necessarily 
banded, and the corresponding covering designs of length 1 higher are unique. At length 4 there 
are four nonisomorphic minimal asymmetric covers, as shown in Table [|, two banded and two 
non-banded. 

Proof. By direct enumeration. The details are omitted. ■ 

At length 7 the unique optimal (and banded) asymmetric covering can be found by deleting 
any coordinate from the following set of 31 vectors of length 8: l 8 ; {0 2 1 2 }1 4 (6); 1 4 {0 2 1 2 }(6); the 
14 vectors of the Steiner system 5(3,4,8); 1 2 6 , 2 1 2 4 , 4 1 2 2 , 6 1 2 . 

4 Or puncturing, cf. p. 28. 
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Table 3: The four nonisomorphic minimal asymmetric coverings 2^(4, 1), of size D(4, 1) = 6. (a) 
and (b) are banded, (c) and (d) are not. 
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Remark. The continuous linear programming problem for the banded case is easily solved, and 
has size exactly 2 n+1 /(n + 2), which is asymptotically 

2 n+1 (--- 2 +0(±)) , (12) 

just slightly worse than ([l0|), 

5. New values for covering numbers 

Let N(n, k, k — 1, M) denote the number of nonisomorphic minimal covering designs C(n, k, k — 1) 
of size M, where of course M > C(n, k,k — 1). The main results of this section are shown in Tables 
| and H 

Table 4: Covering numbers C(n, k, k — 1) for n < 12. Starred entries are new; see text for missing 
entries. 



n/k 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


2 


1 




















3 


2 


1 


















4 


2 


3 


1 
















5 


3 


4 


4 


1 














6 


3 


6 


6 


5 


1 












7 


4 


7 


12 


9 


6 


1 










8 


4 


11 


14 


20 


12 


7 


1 








9 


5 


12 


25 


30 


30 


16 


8 


1 






10 


5 


17 


30 


51* 


50 


45 


20 


9 


1 




11 


6 


19 


47 


66 


a 


84* 


63 


25 


10 


1 


12 


6 


24 


57 


113 


132 





126* 


84 


30 


11 


13 


7 


26 


78 


7 


245 


6 


e 


185* 


112 


36 



Table [| gives the values of C(n, k, k — 1) for n < 12. Starred entries are new, and we have also 
shown that 

(7(14,10,9) = 259 . 

In every case the coverings achieving these bounds were already known, see [|ll]] for references. Our 
contribution has been to show that no smaller covering exist. The five remaining gaps in Table | 
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are at C(ll,6,5) (where 96 < a < 100), C(12,7,6) (165 < p < 176), C(13,5,4) (149 < 7 < 157), 
(7(13,7,6) (257 < 5 < 264), (7(13,8,7) (269 < e < 297), the lower bounds being new, except for 7. 

Table § gives values of N(n, k,k — l, M) together with a brief indication of how they were found. 
An entry such as 

11,8,7: 63(40); 64(1193) (based on [10, 7, 6, 46] ) 

indicates that there are 40 nonisomorphic minimal coverings C(ll,8,7) of size 63, 1193 of size 
64, and that the latter enumeration was based on examining all possible ways to extend minimal 
covering designs C(10, 7, 6) of size < 46. (For any C(ll, 8, 7) of size 63 must puncture to a covering 
C(10, 7,6) which contains a minimal C(10, 7,6) of size < 46.) The symbol ■ at the end of a line in 
the tables indicates that the enumeration of minimal covering designs C(n, k, k — 1) for these values 
of n and k is complete. 

It is worth drawing attention to the gaps that occur just above the parameters corresponding 
to the Steiner coverings C(8,4,3), C(10,4,3), C(ll,5,4), C(12,6,5). For example a C(12,6,5) that 
does not contain the 132-block Witt design must contain at least 137 elements. 

We know of no earlier table of this type, although isolated values have been published. For 
example de Caen et al. @ showed that iV(9, 5, 4, 30) = 3 and iV(10, 6, 5, 50) = 1. The Steiner 



triple systems S(2,3,n) have been enumerated for n < 19 [13]: this gives the number of optimal 
C(n, 3, 2)'s for n = 1 or 3 ( mod 6). Also Steiner quadruple systems S(3, 4, n) have been enumerated 
for n < 15 (see the survey article ]l2|); this gives the number of optimal C(n, 4, 3) for n = 2 or 
4 (mod 6). 

To test isomorphism we generally used the isomorphism subroutines in the Magma computer 
algebra system Q, @, §. 

To compute the entries N(n,k,k — 1,M) in Table || and to establish the new lower bounds 
implicit in Table || we made use of two different branch-and-bound procedures. 

The first procedure branched by selecting one of the uncovered {k — l)-subsets which had the 
fewest remaining /c-subsets which could cover it. Let {x±, X2, ■ ■ ■ , scj} denote the fc-subsets which 
could cover it. The procedure recursively considered the I alternatives {{xj in the covering, xj 
not in the covering, for 1 < j < i — 1}, 1 < i < I}. This branching continued until either every 
{k — l)-subset was covered, or until the lower bound 



E 



j, max^/c ([{>;" G T : v" C v}\) 



(13) 



on the number of available A:-subsets needed to cover the set T of uncovered {k — l)-subsets showed 
that no covering of size M could be obtained from the current branch. 

The second procedure used solutions of the continuous LP problem (see Section 3) to guide 
it. It branched on the fe-subset x whose corresponding variable was closest to ^, considering the 
alternatives "x in the covering" and u x not in the covering". This branching continued until one 
of the following obtained: 

(a) every /c-subset had been placed in or excluded from the covering, 

(b) some {k — l)-subset could no longer be covered, or 

(c) applying the following lemma, where U is the set of available A:-subsets and T is the set of 
uncovered (k — l)-subsets, showed that no covering of size M could be obtained from the current 
branch. The optimal choice of w in this bound is given by the solution to the dual of the linear 
programming relaxation; the program used a discrete, exact approximation to these dual variables 
for its bound. 
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Lemma 1. Let T and U be finite sets equipped with a relation x C y for x € T ', y £U , and let 



w : T 



be a function satisfying 

w(x) < 1 for all y £ U . 

x<ET:x<Zy 

Then any C C U covering T satisfies the lower bound 



(14) 



\C\> 



^2w(x) 



Proof. Let C C U cover T. Then 



\C\ 



> 



> 



yGC x£T:xCy 

^2w(x) 



x£T 



where the first inequality is from (|14|) , and the second is because C covers T. Since |C| is an integer, 
the result follows. ■ 



Remark. The bound ( |i~3D is the special case of the Lemma in which U is the set of all available 
/c-subsets, T is the set of uncovered (A; — l)-subsets, and 



W[V) 



max v:v > Cv (\{v" G T : v" C v}\) 



Because the second program uses a stronger bound, it searched smaller branch- an d-bound trees, 
but since it solved the linear programming relaxation at each node, it took more time per node. 
As a result, the first program was more efficient for "easy" problems, and the second program for 
"difficult" problems (roughly, those in Table 4 in the region bounded by n > 9 and 4 < k < n — 4) . 
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Table 5: Values of N(n,k,k — 1,M), the number of nonisomorphic minimal covering designs 
C(n, k, k — 1) of size M. See text for further details. 

4,2,1: 2(1), 3(1) . 

5,3,2: 4(1), 5(1), 6(1) . 

6,4,3: 6(1), 7(1), 8(1), 10(1) . 

7,5,4: 9(1), 11(2), 12(2), 15(1) . 

8,6,5: 12(1), 13(1), 15(2), 16(3), 17(2), 21(1) ■ 

9,7,6: 16(1), 18(1), 19(2), 20(3), 21(2), 22(3), 23(3), 28(1) ■ 

10, 8, 7 : 20(1), 21(1), 24(4), 25(2), 26(7), 27(5), 28(4), 29(2), 

30(4), 36(1) ■ 

11,9,8 : 25(1), 27(1), 29(2) (based on 10,8,7,21) 

12, 10, 9 : 30(1), 31(1), 34(1), 35(4)(based on 11,9,8,29) 

13,11,10 : 36(1), 38(1), 40(0) (based on 12,10,9,31) 

5,2,1: 3(1), 4(1) . 

6,3,2: 6(1), 7(5), 8(2), 10(1) . 

7,4,3: 12(4), 13(57), 14(139), 15(24), 16(6), 17(1), 20(1) ■ 

8, 5, 4 : 20(6), 21(263), 22(7340)(based on 7,4,3,13) 

9, 6, 5 : 30(2), 31(16), 32(863)(based on 8,5,4,21) 

10, 7, 6 : 45(20), 46(609)(based on 9,6,5,32) 

11, 8, 7 : 63(40), 64(1193) (based on 10,7,6,46) 

12, 9, 8 : 84(4), 85(46), 86(1423) (based on 11,8,7,64) 
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Table 5: (continued) 



6,2,1: 3(1), 4(2), 5(1) . 

7,3,2; 7(1), 9(14), 10(40), 11(60), 12(7), 13(1), 15(1) . 

8,4,3: 14(1), 17(13) 

9, 5, 4 : 30(3), 31(18), 32(459)(based on 8,4,3,17) 

10, 6, 5 : 50(1), 52(4), 53(56), 54(880)(based on 9,5,4,32) 

11, 7, 6 : 84(3), 85(0)(based on 10,6,5,54) 

12. 8, 7 : 126(3), 127(2), 128(0)(based on 11,7,6,84) 

13. 9, 8 : 185(1), 186(0)(based on 12,8,7,127) 
14, 10, 9 : 259(l)(based on 13,9,8,185) 

7,2,1: 4(1), 5(2), 6(1) . 

8, 3, 2 : 11(5), 12(145)(based on 7,2,1,4) 

9, 4, 3 : 25(77), 26(5562), 27(538969) (based on 8,3,2,12) 

10, 5, 4 : 51(40), 52(3354) (based on 9,4,3,26) 

8,2,1: 4(1), 5(2), 6(3), 7(1) . 

9, 3, 2 : 12(1), 13(1), 14(64)(based on 8,2,1,4) 
10, 4, 3 : 30(1), 33(43)(based on 9,3,2,13) 

11, 5, 4 : 66(1), 70(78)(based on 10,4,3,30) 
12, 6, 5 : 132(1), 137(87)(based on 11,5,4,66) 

9,2,1: 5(1), 6(3), 7(3), 8(1) . 

10. 3, 2 : 17(58)(based on 9,2,1,5) 

11.4, 3 : 47(> 95970)(based on 10,3,2,17) 
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